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Abstract. We present photometry and spectroscopy of the afterglow of GRB 030323. VLT spectra of the afterglow 
show damped Lya (DLA) absorption and low- and high-ionization lines at a redshift 2=3. 3718±0. 0005. The 
inferred neutral hydrogen column density, log N{Hl)=21.90±0.07, is larger than any (GRB- or QSO-) DLA Hi 
column density inferred directly from Lya in absorption. From the afterglow photometry, we derive a conservative 
upper limit to the host-galaxy extinction: Av<0.5 mag. The iron abundance is [Fe/H]=-1.47±0.11, while the 
metallicity of the gas as measured from sulphur is [S/H]=-1.26±0.20. We derive an upper limit on the H2 
molecular fraction of 2Ar(H2)/(2iV(H2)-|-Af(Hl))^10"'^. In the Lya trough, a Lya emission line is detected, which 
corresponds to a star- formation rate (not corrected for dust extinction) of roughly 1 M© yr~^. All these results 
are consistent with the host galaxy of GRB 030323 consisting of a low metallicity gas with a low dust content. We 
detect fine-structure lines of silicon. Sin*, which have never been clearly detected in QSO-DLAs; this suggests 
that these lines are produced in the vicinity of the GRB explosion site. Under the assumption that these fine- 
structure levels are populated by particle collisions, we estimate the Hi volume density to be uhi = 10^ — 10* 
cm~"^. HST/ACS imaging 4 months after the burst shows an extended AB(F606W)=28.0±0.3 mag object at a 
distance of 0'.'14 (Ikpc) from the early afterglow location, which presumably is the host galaxy of GRB 030323. 
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1. Introduction 



Table 1. Log of UT4/FORS2 spectroscopic observations 



Damped Lya (DLA) absorbers, conventionally detected 
in Quasi-Stellar Object (QSO) spectra, are absorption- 
line systems that have a column density of N(Hi) > 
2x10^° atoms cm~^, as determine d from the damp- 
ing wings of the Lya line (e.g. I Wolfe et al.1 llQSfit 
iTurnshek et al.lll989() . DLA systems are believed to con- 
tain the bulk of the neutral hydrogen at high red- 
shift and to be the gas reservoir from wh ich the stars 
at the present epo ch are produced (e.g. IWolfel Il987t 
iLanzetta et aLlHoQlj) . Numerous high-resolution spectro- 
scopic studies have extracted detailed info rmation about 
the metallicity (e.g.|P rochaska et al."2003a V the kinemat- 
ics l|Prochaska fc WoTfgilQOTtlLedoux et al.il998t). and the 



dust and H9 contents l)Petitiean et al.ll2n0f][ iLedoux et al.l 

1200.1(1 of DLA systems as a function of redshift. Despite 
intensive searches, only a han dful of DLA counte rparts 
have been detected so far (see iMoller et al.ll2002bl) : link- 
ing DLA systems with galaxy types has therefore proven 
difficult: some ad vocate large, disk-forming galaxies (e.g. 
IWolfe et ahllTflfli IProchaska Wolfelll 997ll. others su 
gest they are faint, gas- rich dwarfs l)Haehnelt et alJll99i 

Gamma-ray burst (GRB) afterglows are, just as 
QSOs, bright and distant sources. For instance, the 
spectacular GR B 990123 was dete cted at the 9*^ vi- 
sual magni tude llAkerlof et al.lll999t) while it was l ocate d 
at z=1.6 l|Kulkarni et all Il999t lAndersen et all Il999|) . 
However, the afterglow brightness in general fades very 
rapidly in time (roughly flux oc time~^). The current 
afterg low redshi fts range from z=0. 169 ijGreiner et al.l 
12003(1 to z=4.5 l) Andersen et alJl2000() . Moreover, GRBs 
are associated with massive-star formation: the dis- 
covery of a su pernova (SN) spectru m similar to that 
of SN1998 bw l)Galama et all Il998bl) superimposed on 
the GRB 030329 aft erglow spectrum l(Stanek et al.ll2003t 
iHiorth et alJ l2003bll provided strong evidence that at 



least some of the long-duration (^ 2s) G RBs are caused 
by the core collapse of ma ssive stars l|Wooslevl Il993t 
iMac.Fadven Woosleviri99^. 

The discovery of a damped Lya (DLA) absorption line 
at the burst redshift in the spectra of several GRB after- 
glows il.Ten sen et aLlEnml iFvnbo et aljEoml iHiorth et all 
l2003a|l is consistent with the massive-star progenitor sce- 
nario: they indicate a high neutral hydrogen column den- 
sity origin in the host galaxy, presumably a star-forming 
region. However, the signal-to-noise ratio at the location 
of the DLA absorption line in the spectra is fairly low 
in these cases, much lower than for typical QSO-DLAs. 
We here present afterglow spectra of the high-redshift 
GRB 030323, which unambiguously demonstrate a GRB- 
DLA, with a column density exceeding that of any (QSO- 
or GRB-) DLA measured so far using Lya in absorption. 



Send offprint requests to: pvreeswi@eso.org 

* Based on observations collected at the European Southern 
Observatory, Chile, by GRACE (Gamma-Ray Burst Afterglow 
Collaboration at ESQ), under program ID 70.D-0523(B), and 
with NASA's Hubble Space Telescope. 
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These spectra allow us to measure the metallicity of the 
host environment and obtain an upper limit on the molec- 
ular fraction, i.e. measurements that are routinely per- 
formed for QSO-DLAs, but that are still unique for GRB 
hosts. Although the GRB-DLA sample is still very small, 
we compare them with QSO-DLAs in two aspects: their 
Hi column density and their metallicity. 

GRB 030323 was detected on 23 March 2003 at 

with a flu- 



21:57 UT by HETE l|Graziani et al 
ence of 1.1x10"^ ergs cm~^ (30-400 keV), and a du- 
ration of 26 seconds. Following the HETE localization, 
the optical counterpart was discovered 7.6 hours af- 
ter the bur st at R.A. ll '^ 06"'0ff :38, Decl. -21°46'13'.'3 
(J2000) by iGilmore et ail l|2003(l . with a brightness of 
R=18.7. Our team reported a preliminary redshift of 
z=3.372 (^yrecswiik et al. 2003), which is current l y the 



third highest redshif t for a GRB l(KulkarnietHI Il998t 
lAndersen et ^Emh. 

This paper is organized as follows: in Sect. [3 we de- 
scribe the data reduction of both the spectroscopic and 
imaging observations. In Sect. |21 we present the light 
curves and infer an upper limit on the rest-frame opti- 
cal extinction. We measure the equivalent widths of the 
absorption lines and determine the burst redshift in Sect. 
21 An Hi column density model is fitted to the damped 
Lya line in Sect.|3 and we analyze the spectra in more de- 
tail in Sect. El to derive the metallicity and an upper limit 
on the molecular hydrogen (H2 ) fraction. The detection of 
Lya in emission is presented in Sect.[71 and we report on 
the detection of the probable host galaxy of GRB 030323 
in HST/ACS imaging data in Sect.|S| In the final section, 
we close with a general discussion of all these results. 

2. Observations and data reduction 

The spectroscopic observations of GRB 030323 were 
performed with the Focal Reducer Low Dispersion 
Spectrograph 2 (F0RS2) at unit 4 (Yepun) of the 
Very Large Telescope (VLT) at the European Southern 
Observatory (ESO) at Paranal, Chile. The imaging obser- 
vations were performed with a variety of telescopes and 
instruments. Tables ^ and [21 show the spectroscopy and 
imaging observation logs. 

The images and 2-D spectra were bias-subtracted and 
flat-fielded in the usual manner, mostly within IRAF^. 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
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Table 2. Log of imaging observations 
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15 


0.7 


F0RS2 


Jul 5.993 


> 25.6 (3a) 


R 
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The magnitudes have not been corrected for Galactic ex- 
tinction, and the errors do not include the uncertainty in 
the absolute calibration (see text). 

Combinations of telescopes and instruments; USNO 
Flagstaff Station 1.0m with 2kx2k Tek CCD, O'.'68/pixel; 
Calar Aho 2.2m and CAFOS with Ikxlk SITe CCD, 
O'.'29/pixel; ESO/Danish 1.54m and DFOSC with 2kx4k 
EEV CCD, binned to O'.'78/pixel; Yepun and F0RS2 with 
two MIT CCDs of 4k X 2k, binned to O'.'25/pixel; Yepun and 
NACO with the Aladdin InSb Ikxlk detector and S54 
camera, O'.'054/pixel; SARA 0.9m and Apogee Ap7 CCD 
camera with 512x512 array, O'.'7/pixel; UKIRT and UFTI 
with a Ikxlk HgCdTe array, O'.'091/pixel; Gemini South 
and acquisition camera with IKxlK EEV CCD, binned to 
O'.'23/pixel; HST/ ACS and WFC detector with two 2kx4k 
SITe CCDs, 0'.'05/pixel. 



Following this, the spectra were cosmic-ray cleane d usin g 
the L.A. Cosmic program written by Van Dokkunj l)200l|) . 
Each spectrum was extracted separately, and wavelength- 
calibrated using an HeNeAr lamp. The error in the wave- 
length solution was of the order of O.lA for the low reso- 
lution 300V and 3001 grisms (with A/ A A of 440 and 660, 
respectively), and O.OSA for the intermediate resolution 
grisms 1400V and 1200R (with A/AA - 2100). 

Flux calibration was performed using the standard 
LTT3864, and the slit losses were determined for each 
grism by fitting a Gaussian along the spatial direction of 
the summed 2-D spectra, every 4 pixels across the entire 
dispersion axis (i.e. summing 4 columns before performing 
the fit). The resulting Gaussian full width at half maxi- 
mum (FWHM) was then compared to the slit width to ob- 
tain the slit loss (i.e. the fraction of the surface underneath 
the Gaussian fit that is outside the slit width) along the 
dispersion axis. The slit loss profile was then fitted with a 
polynomial to correct the spectra. All GRB spectra were 
taken with a 1" slit, and slit losses were as high as 60% for 
the blue part of the 300V grism, while around 15% for the 
other grisms. The standards were observed with 5" slits, 
and therefore do not suffer from slit losses. We note that 
Yepun/F0RS2 contains a linear atmospheric dispersion 
compensator (LADC) in the light path, which minimizes 
any colour-differential slit losses up to a zenith distance of 
45°. However, the last three 1200R spectra were taken at 
an airmass above this limit, and are therefore affected. 

Finally, the first night's 300V+300I and second night's 
1400V+1200R spectra were combined into two full spec- 
tra and corrected for Galactic extinction ijSchlegel et alJ 
1998). The scahng between the blue and red parts in the 
combining of the 300V with 3001 spectra, and 1400V with 
1200R, was determined from an overlapping region free 
of lines, and amounted to 3.5% of the continuum level be- 
tween 300V and 3001, and 2% for 1400V and 1200R, which 
can be explained by the fading of the afterglow between 
the epochs at which the spectra were taken. 

Comparison between the spectroscopy and imaging 
absolute calibration shows that the spectroscopic flux is 
roughly 30% below the BVRI photometric measurements 
(i.e. not colour dependent) that were taken around the 
same epoch. We have not corrected the spectrum in Fig. 
Ofor this difference. 



3. Photometry 

The magnitudes of the afterglow (see Table were deter- 
mined through aperture photometry, relative to a set of 
reference stars. The aperture radius was set to the average 
FWHM of the point spread function of stars on the image 
(column 5 in Table[2Jl. The early July upper limits are 3(7, 
and also for an aperture radius the size of the seeing disk. 
The reference stars in t he optica l were tied to the absolute 
calibration provided bv lHendenI l|2003^. This was found to 
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1 10 100 

Days since burst (March 23.915, 2003 UT) 



Fig. 1. BVRI JHK light curves of GRB 030323. The magnitudes have been corrected for the Galactic extinction of 
Eb-v=0.049 ijSchlegel et alJlTgOSl) : note that in Table |21 they have not been corrected. The dashed lines are simple 
power law fits to the observations before day 2, and extrapolated to later epochs. 



be consistent with our own calibration on two different 
nights, but during which only one standard field was ob- 
served. For the near-infrared filters we used the calibration 
provided by the 2MASS^, as none of the nights on which 
the near-infrared data were taken were photometric. The 
magnitude errors listed in Tableware a combination of the 
Poisson error of the afterglow measurement and the scat- 
ter in the reference star magnitudes; they do not include 
the error in the absolute calibration. For the optical band- 
passes, we determined this absolute error by calculating 
the average differen ce of two calibr ations on two different 
nights provided by 'HenderJ l)2003l) for several stars; we 
find the following magnitude errors: 0.02 (B), 0.03 (V), 
0.02 (R), and 0.f2 (1). For the infr ared filters th e abso- 
lute calibration error is provided bv lCohen et alJ ((20031) . 
which amount to 0.02 mag for J, H and K. No correc- 
tion has been attempted for the fact that the observations 
were taken with filters from different systems. The mag- 
nitudes listed in Table |2 have not yet been corrected for 
the Galactic foreground absorption. 

Fig. n shows the light curves of GRB 030323 in sev- 
eral wavelengths. The light curves can be fit neither with 
a single power law, nor with a broken power law. This is 
clear from performing simple power law fits to the early 
epochs in the optical bands, and extrapolating these to 
later times. These fits are shown by the dashed lines in 

^ see |http: / /www.ipac.caltech.edu/2ma8s/| 



Fig.n The V band light curve shows the clearest deviation 
in between day 2 and 3 after the burst. However, after day 
3 the observations are fairly well described again by the 
extrapolation of the early slope, even the R band measure- 
ment around day f 0. The V band extrapolation underes- 
timates the late-time HST point by roughly a magnitude, 
which suggests that the afterglow has become fainter than 
its host galaxy at this epoch. 

Several afterglows have displayed deviations 
from the commo n smooth power law decay, such as 
GRB 970 508 (e.g. 'Galama et al.' 'f 998a>), GRB 000301C 
(e.g.jMasetti et all 2000). GRB 02f 00 4 (e.g.lHolland et all 
20031). and G RB 030329 (e.g. IPrice et alJ l2003^ . 
Garnavich et al.l |2000) have suggested that a mi- 
crolensing event caused the deviations in GRB00030fC. 
In the context of the fireball model, the deviations in 
the case of GRB 021004 are interpreted as due to a 
variable external density l|Lazzati et al.l l2002bj) . and for 
GRB 030329 they are interpre ted as due to refr eshed 
shocks from the inner engine llGranot et al.ll2003^ . We 
note that the cannonball model off ers an alternative 
explanat ion for these observations ijPado et alJ I2002L 
l2003alb|) . For GRB 030323, we only study the global 
properties of the light curve, and compare them to the 
light curve decay and spec tral slope values as pre dicted 
by the fireball model (e.g. ISari et all Il99l Il999|) . with 
the aim of constraining the host-galaxy extinction. 
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As the late-time afterglow behaviour is not clear, we 
only use the simple power law fits to the early optical 
data. The inferred optical temporal decay indices are con- 
sistent with one another, with an average decay of aopt=~ 
1.56±0.03 (using the convention: Fj^(t) (xt"i^^). The near- 
infrared slopes have similar values as the optical ones: aj 
= -1.81±0.39, an = -1.42±0.14 and = -1.46±0.31. 
However, these may be affected by the "bump" in between 
day 2 and 3, if it is achromatic. 

At several epochs after the burst, observations in at 
least two filters were performed around the same epoch. 
This allows us to construct broad-band spectral energy 
distributions (SEDs) and fit them to obtain the opti- 
cal to near-infrared spectral slopes at these epochs. Note 
that we discard the B and V bands, as these are attenu- 
ated by the Lya line and forest absorption. For day af- 
ter burst 0.40, 1.18, 2.17, 2.44, 3.27, and 4.50, we ob- 
tain: /3(0.40)= -1.05±0.64 {Xrcd=^)^ /?(1.18)= -0.83±0.61 
(X?ed=0), /3(2.17)=~1.09±0.07 (x?,d=4-5), /3(2.44)=~ 
0.80±0.05 (x?^d=l-7), /3(3.27)=-1.2±1.0 (xld^O), and 
/?(4.50)=-0.82±0.12 (x?ed=4-5). For the epochs with 
Xred~^' observations in only two filters are available. 
Except for the fit value at day 2.17, these values are con- 
sistent with being constant, and the weighted mean and 
its error is /3obs= -0.89±0.04. 

We now compare these observed spectral slopes with 
the ones predicted by the fireball model to obtain an es- 
timate of the host-galaxy extinction. An important as- 
sumption that we make in estimating the optical extinc- 
tion is that the intrinsic afterglow spectrum is a power 
law, which is a consequence of the fireball theory for GRB 
afterglows. The fireball theory has be en quite successfu l 
in explaining the observations (but see lDado et alJl2002|) . 
The predicted spectral slope depends on the assumed cir- 
cumburst densit y profile being either c onstant or that of a 
stellar wind fsee lLi fc Chevalie3l200lh . whether th e light 
curve is in the jet regime or not ( see LSari et aljfl999il . 
and whether the cooling break (see ISari et alJll998(l has 
already passed the optical wavebands (cooling regime) or 
not. Considering all these possibilities, the predicted spec- 
tral slope ranges from /3cxp=(a+l)/2=-0.28±0.01 (wind 
or constant density medium, post jet-break and non- 
cooling regime) to /3cxp=(2a-l)/3=-1.37±0.01 (wind or 
constant density, pre jet-break and cooling regime), with 
a=-1.56±0.03. We have assumed that these relations be- 
tween a and P are also valid for a power la w index of the 
electr on energy distribution, p < 2 (but see lPai fc Chend 

l2nnij^ . 

All intrinsic spectral slopes shallower than the ob- 
served slope of -0.89 =b0.04 leave some room f or host- 
galaxy extinction (see 'Ram aprakash et al.lll99i ). as any 
host-galaxy extinction results in a steepening of the in- 
trinsic slope. We conservatively take -0.28±0.01 to be 
the actual slope, to obtain an upper limit on the h ost- 
galaxy extinction. Using the extinction-curve fits of IPeil 
l)l992t) for the Milky Way (MW), and the Large- and 
Small Magellanic Clouds (LMC and SMC), we itera- 
tively find Ay that fits best with the expected spec- 



K H ^ J ^ I ^ R 

MW 




Fig. 2. Milky Way (MW), Large- and Small Magellanic 
Cloud (LMC and SMC) extinction-curve fits to the spec- 
tral energy distribution (SED) of the afterglow at day 
2.44 after the burst, assuming that the intrinsic SED is 
a power law with slope /3= -0.28 (see text). These fits 
result in the conservative upper limits to the host-galaxy 
extinction: Av(MW)<0.50 mag, Av(LMC)<0.25 mag and 
Av(SMC)<0.16 mag. 

tral slope of -0.28 (held fixed in the fit) for the en- 
ergy distribution at 2.44 days after the burst. We find 
Av(MW)=0.50 mag (^2^^=1.5), Av(LMC)=0.25 mag 
(x?cd=l-8) and Av(SMC)=0.16 mag (x?cd=2.3). These 
fits are shown in Fig. [21 We note that if we would 
have assumed that the light-curve break occurred after 
day 1.4, which is hkely, then the expected spectral slope 
would be /?cxp=(2a-l-l)/3=-0.71, and the extinction val- 
ues would decrease to Av(MW)=0.09 mag (^^^^=1.3), 
Av(LMC)=0.04 mag (x?ed=l-4) and Av(SMC)=0.02 mag 
(Xrcd— When we do not fix the intrinsic slope at a par- 
ticular value (but still assume that the intrinsic spectrum 
is a power law), we find spectral slopes ranging from -0.80 
to -0.67, and Ay from to 0.12 mag. Therefore, a very 
conservative upper limit on the host galaxy extinction is: 
Av<0.50 mag. 

4. Absorption-line measurements and redshift 

The combined 1400V-K1200R spectrum of GRB 030323, 
corrected for Galactic extinction, is shown in Fig. |3| The 
most obvious feature is the very broad absorption fine 
around 5300A, which can be identified as Lya. Redward 
of Lya several metal absorption lines are detected, and 
to the blue the intervening Lya forest is present. We also 
show the 1(7 Poisson error spectrum. 

After normalization of the spectrum with a high-order 
(25) polynomial, we measured all possible lines with splot 
in IRAF, summing the equivalent width (EW) of the indi- 
vidual pixels, and determined the line center. In case two 
lines were clearly blended, we used splot to deblend, using 
Gaussian line shapes, and forcing a single FWHM for both 
lines. The EWs and centers of the lines above 5a signifi- 




Fig. 3. Combined 1400V+1200R spectrum of GRB 030323, including the Poisson error spectrum. The inset on the left shows the normalized spectrum with the 
hydrogen column fit to the damped Lya line, including the la errors. This column density is currently the highest for any DLA system measured using Lya in 
absorption (see Fig.0J. Redward of Lya we detect numerous metal-absorption lines (see Tabled), whose average redshift (in this high-resolution spectrum) is 
2=3.3716 ± 0.0005. In the inset on the right, we zoom in on the CiV doublet, which is split into separate components. In the Lya trough, Lya is detected in 
emission. Lines that could not be identified, atmospheric absorption lines (indicated with an ©), and lines belonging to the tentative absorber at z—lAl, are 
indicated with a dashed line instead of a solid one. 
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cance are tabulated in Tabled for both the low- (Ir) and 
high-resolution (hr) spectrum, along with the error in the 
EW, the line identification and the line redshift. These line 
redshifts are determined from the high-resolution spec- 
trum, unless the line is not covered by this spectrum (i.e. 
above 73f OA and below 4560A - see Table^l. The error in 
the EW is determined by: AEW = AA v^Sa?, where the 
error spectrum (a) is normahzed by the high-order polyno- 
mial continuum fit of the object spectrum, and is summed 
over the same pixels (i) that were used to measure the line 
EW. A A is the number of A per pixel, which is 3.2 and 
0.64 A/pixel for the combined 300V+I and 1400V-H1200R 
spectra, respectively. 

Using all lines that are detected above 5a significance 
and that could be identified, we find z(lr)=3.3728±0.00f f 
and z(hr)=3.3716±0.0005, for the low- (h) and high- 
resolution (hr) spectrum, respectively. We adopt the 
weighted mean of these values as the redshift of 
GRB 030323: z=3.3718±0.0005. There is no doubt that 
this is the redshift of GRB 030323 and not that of a 
chance foreground galaxy, since otherwise Lya forest lines 
redward of the DLA line would have been detected. We 
find that the strongest lines in the red part of the high- 
resolution spectrum, CivAA1548, 1550, are split into sep- 
arate components (see the inset of Fig. with a velocity 
difference of 130±60 km s^^ (where the error on the wave- 
length determination of the lines in the blend is 1 A) . Such 
a velocity spread is consistent with the absorption taking 
place in separate regions in the host galaxy. 

Several lines that are detected above 5a can not be 
identified, some of which correspond to significant lines 
in the standard star spectra. Most of these are imprinted 
on the spectrum by the Earth's atmosphere. Several lines 
not belonging to the GRB host galaxy can in principle 
be identified with Fell AA2344.2, 2600.1 (although an at- 
mospheric line is detected at the latter wavelength) and 
Mgii AA2796.3, 2803.5, all around z=1.41. However, the 
oscillator strength of Fe II A2382 is three times larger than 
that of the detected Feii A2344, but this line is not de- 
tected. Moreover, both Fell lines are s tronger than Mgii, 
which is usually not observ ed (e.g. iBoisse & BergeronI 
ll985HSteidel fc Sar genii 992l) . Therefore, we consider the 
existence of this foreground absorption system to be highly 
uncertain. 

The optical/UV flash of the GRB is expected to alter 
its immediate environment, possibly leading to a change 
in absorption-line strengths as a function of time (see 
Perna & Loeb' 'l998^ IVreeswiik ei~an l200l[ IPerna et alJ 
2003) . Comparing the low- and high-resolution equivalent 



Table 3. Lines detected above 
high-resolution (hr) spectra. 



5a in the low- (Ir) and 



A(lr) EW,ost(lr) A(hr) EWrcst(hr) 



widths in Table 13 shows that none of the lines detected in 
both spectra are significantly varying {3a). 

5. H I column density 

We have fitted a power law continuum to the high- 
resolution spectrum over the wavelength range 5870- 
7000 A and determined a power law slope (in F^) of - 
0.76±0.09 (fitting only the first 1200R spectrum, as the 



4484 
5317 



5483 
5511 



9(1)^ 



5692 
5706 

5706 
5727 
5838 

6095 
6135 
6280 

6676 
6706 



6770 
6782 
6875 
6900 
7034 



7307 
7907 
8112 
8145 
8237 
8802 
8859 
8970 
9019 



0.44(6) 
1.07(5) 



5530 0.76(5) 



5637 0.33(5) 



0.75(6) 
1.08(6) 

1.08(6) 
0.26(5) 
1.71(7) 

1.19(6) 
0.98(6) 
0.51(3) 

0.69(4) 
0.54(4) 



1.42(4) 
1.28(4) 



0.51(4) 



1.08(5) 
0.58(5) 
0.45(5) 
0.23(4) 
0.34(4) 
0.55(6) 
0.82(9) 
0.73(11) 
0.77(5) 



5316.9 
5415.8 
5433.5 
5465.4 
5467.4 
5481.1 
5505.0 
5509.7 

5529.2 
5591.5 
5595.7 
5624.4 
5634.7 
5649.2 
5652.0 
5692.5 
5702.4 

5709.7 
5723.9 
5833.7 
5838.9 
6092.9 
6132.4 
6278.0 
6629.2 
6673.8 
6703.7 
6743.0 
6754.2 
6760.6 
6766.9 
6769.9 
6778.2 
6781.0 
6870.4 
6876.5 
6891.9 
6902.4 
7031.7 
7043.6 
7082.0 
7160.1 
7202.7 
7303.5 



0.31(3 
0.20(2 
0.15(2 
0.21(3 
0.29(2 
0.33(2 
1.00(2 

0.75(2 
0.11(2 
0.10(2 
0.11(2 
0.16(2 
0.16(1 
0.15(1 
0.69(2 
0.70(2 

0.32(2 
0.29(2 
0.84(3 
0.79(2 
1.08(4 
1.15(4 
0.42(3 
0.14(2 
0.75(3 
0.34(3 
0.14(2 
0.13(2 
0.14(3 
0.99(3 
0.51(3 
0.84(3 
0.53(3 



0.48(3) 
0.17(3) 
0.15(3) 
0.20(3) 
0.17(3) 
0.99(7) 



Ly/3A1025.7 3.3715 

Lyaem. A1215.6 3.3736 

NvA1238.8 3.3717 

NvA1242.8 3.3720 

SIIA1250.5 3.3719 

S11AI253.8 3.3716 

S11AI259.5 3.3707 

SiiiA1260.4 3.3713 

orFellA1260.5 3.3709 

Siii*A1264.7 3.3718 

©,CiA1280.1? 3.3679 



FellA2344.2? 1.4110 

O1AI302.I 3.3716 

SiiiA1304.3 3.3718 

andOi* A1304.8? 3.3702 

Oi**A1306.0? 3.3718 

Siii*A1309.2 3.3718 

C11AI334.5 3.3713 

Cii*A1335.7 3.3714 

SiivA1393.7 3.3716 

SiivA1402.7 3.3716 

®,FellA2600.1? 1.4145 

SiiiA1526.7 3.3714 

Siii*A1533.4 3.3717 

MgiiA2796.3? 1.4114 

MgiiA2803.5? 1.4092 

CivA1548.2 3.3708 

CivA1548.2 3.3728 

CIVA1550.7 3.3708 

CIVA1550.7 3.3727 



FeiiA1608.4 3.3717 
FeiiA1611.2 3.3716 



A1IIA1670.7 3.3713 

SiiiA1808.0 3.3734 

AI111AI854.7 3.3738 

AI111AI862.7 3.3724 



ZniiA2026.1 3.3725 

e 

3,ZniiA2062.6 3.3726 



Ly/3 is located in the Lya forest, which causes it to be 
blended with forest lines, and the continuum placement is 
highly uncertain, resulting in a large EW error. 
The lines at AA5652.0, 6278.0, 6743.0, and 6754.2 can in 
principle be identified with Fell and Mgll at Z=1A1; how- 
ever, strong absorption is then expected from Fell A2382, 
which is not detected. 

The lines for which significant absorption was also detected 
in the standard star spectrum are marked with a ®. 
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other 1200R spectra suffer from colour-dependent slit 
losses). This value is in agreement with the red slope of 
the 3001 spectrum and with the slope of the photom- 
etry measurements (see Sect. We used an extrapo- 
lation of this power law to blue wavelengths in order 
to normalize the entire spectrum. The resulting average 
flux decrement in the Lya forest bet ween Ly/3 and Lya: 
Da = 1 — (observed) /Fj, (intrinsic) ijOke fc Korvcanskvl 
Il982(l . that we obtain is Da=O.U ±0.04. This decrement, 
which is due to intervening hydrogen systems, is consistent 
with that obs erved in QSO lines of sight at the redshift of 
GRB 030323 llCristiani et alJll993l) . 

A fit to the strong Lya absorption (using the DIPSO 
package within Starhnk) yields log N(H i)=21.90±0.07. 
This fit is shown in the inset of Fig. O This is the 7th 
GRB for which a neutral hydrogen column has been de- 
termined from the afterglow spectrum, and GRB 030323 
happens to have the highest H i column density measured 
so far. Fig.^shows a comparison of the H i column density 



Curran et al."2002') and GRB-DLAs ( 


Jensen et al.ll200ll 


Fvnbo et al...2001i- iHiorth et al.ll2003a 


, Jakobsson et al. 



2004, in prep., and this paper). For completeness, we also 
show the two GRBs for which Lya was detected but which 
do not qualify as a DLA system: GRB 0112 11 (Vreeswiik 
et al 2004, in prep.) and GRB 021004 l|Moller et all 
l2002a^ . It is quite striking that out of 7 GRB afterglows 
for which Lya was red-shifted into the observable spec- 
trum, 5 show evidence for a high column density DLA 
system. This clearly demonstrates that GRBs explode in 
either galaxies, or regions within galaxies with high ncu- 




Fig. 4. Histogram of the column densities of DLA sys- 
tems measured through the damping wings of Lya dis- 
covered in the spe ctrum of a backgro und QSO (compi- 
lation taken from ICurran et aLll2002j) . The shaded his- 
togram shows measurements in GRBs for which the red- 
shift was large enough to detect Lya. Out of 7 GRBs, 5 
show neutral hydrogen column densities above the DLA 
definition of 2x10^° atoms cm'^ (log N(Hi)=20.3). The 
host of GRB 030323 contains a column density larger than 
in any observed (GRB- or QSO-) DLA system. 



tral hydrogen column densities. The Hi gas responsible 
for these large columns could be related to the site of the 
GRB explosion, e.g. part of the massive-star forming re- 
gion in which the GRB occurred, but could also be gas 
that is not associated with the GRB, further away in the 
host gala xy. We performed a Kolgomorov-Smirnov (KS) 
test (e.g. IPress et 311119921) to estimate that the proba- 
bility that both samples are drawn from the same par- 
ent distribution is 0.0006. Moreover, in this comparison 
with QSO-DLAs, the GRB-DLA Hi column densities are 
in fact lower limits as the GRB itself occurs within the 
galaxy that is associated with the DLA system; if the GRB 
sightlines would have been probed with background QSOs, 
their column densities would have been on average a factor 
of two larger, which would shift the GRB column densities 
in Fig.01by 0.3 dex upward. However, GRB 011211 would 
then move into the GRB-DLA sample resulting in a only 
a slight decrease in the above-mentioned KS probability. 

6. Metallicity and H2 content 

Although there are many metal lines observed in the spec- 
trum of this GRB, most of them are saturated in the in- 
termediate resolution spectrum. We have identified only 
2 sets of lines as potentially unsaturated, based on their 
small (< O.4A) rest-frame equivalent widths (EWs). These 
are the Sii A1250, 1253, 1259 triplet and Fell A1611. Two 
of the Sii lines (A = 1250,1259) show signs of blend- 
ing, evidenced by a weak component that broadens the 
A1250 line in its blue wing, and a strong interloper red- 
ward of the A1259 line. Sii A1253 appears as an unre- 
solved single component. We measured the observed EWs 
of Sii A1253 and Feii A1611 to be 1.25A and 0.72A 
(EWrcst= O.29A and 0.16A), respectively. In the opti- 
cally thin limit, these correspond to the column densities: 
log N(Sii)=15.3 and log N(Fe ii) = 15.7, and abundances 
[S/H]=-1.8 and [Fe/H]=-1.7. In this conversion from col- 
umn dfinsit^_to_abundanc^^ assumed the Solar values 
from^ Grevesse fc Sauvall l)l998(l . and no correction for ion- 
ization; i.e. we assumed that the column densities of S 11 
and Fell are equal to the total column densities of S and 
Fe, as the singly ionized state of both of these elements 
should be the dominant one in a region with such a high 
H I column dens ity. This has been motivated the oretically 
for QSO-DLAs ('ViegaE lll995tlVladilo et ^ ''2001'): we here 
assume that ionization corrections are also negligible in 
GRB-DLAs. 

In a second step, we performed a simultaneous one- 
component fit to the lines Sii A1250, 1253, 1259 and 
Fell A1611, taking into account line blending and a 
range of broadening parameters. Fig. shows the re- 
sulting fits. We find log N(S ii)=15.84±0.19 and log 
N(Feii) = 15.93±0.08, together with a (turbulent) broad- 
ening parameter b=35±10 km s^^. Comparing these col- 
umn densities with those from the optically thin limit 
approximation shows that the Sii A1253 and Fell A1611 
lines are slightly saturated. These measured column den- 
sities correspond to the abundances: [S/H]=-1.26±0.20 
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Fig. 5. Simultaneous one-component fit to low-ionization 
metal lines with little saturation. 



and [Fe/H]=-1.47 ±0.11. Therefore, although the N(Feii) 
is large (see also iSavaglio et al. 2003), the [Fe/H] is only 
marginally higher than that of QSO-DLAs at this redshift: 
the mean \Fe/B] of 26 QSO-DLAs with 3.0< z <3.5 from 
fProchaska_et_aL ( ^OOSafl is -1.83, with a scatter of 0.35. 

Since essentially all DLA systems observed toward 
QSOs have [S/Fe] greater than zero fsee lLooez fc Ellison! 
I2OO3 . for a recent compilation and discussion), our mea- 
surement of [S/Fe]=0.21±0.23 in GRB 030323 does not 
represent a very stringent constraint on a possible a- 
element overabundance. Moreover, as iron is a known 
dust-depleted element, there could be a small correction to 
this ratio due to dust depletion. Interestingly, a tendency 
toward high values of [Si/F e]'^ has been found in other 
GRBs llSavadio et alJl2003l) . as expected in cases where 
massive-star formation has recently deposited metals into 
the ISM. 

We have examined the GRB 030323 spectra for pres- 
ence of H2 absorption lines, but these are not detected 
(for a list of lines and t heir oscillator strengths, see 
iMorton fc Dinersteinlll976|) . The location of possible H2 
lines at Zabs=3.3716 is actually observed for the L=0 to 3 
Lyman bands of H2. Of these, only the expected location 
of the L=2 band is clear of blending with Lya forest lines. 
Because of the low resolution of the spectra, two ranges of 
possible broadening parameters bng (with bna ^ bmctais) 
were considered to perform trials of Voigt-profile fitting 
of both the J=0 and 1 hues (namely: H2 L2-0 R(0), L2-0 
R(l) and L2-0 P(l)). We find the following upper hm- 
its: (1) for the range 10 km s^^ < hn^ ^ 50 km s~^: 
log N(J=0)<14.5 and log N(J=1)<15.5, and (2) for the 



range 1 km s^^ < bn^ < 10 km s^^: log N(J=0)<14 
and log N(J=1)<18. Therefore, strictly speaking the de- 
rived upper limit on the mean molecular fraction of the 
system (i.e. GRB environment -I- host galaxy) is: f = 
2N(H2)/(2N(H2)-hN(Hi)) < 2xl0i7(2xl0i8-H7.9xl02i) 
= 2.5x10-'* with N(Hi)=7.9xl02i cm'^ (i.e. the 
above case 2). However, under the assumption that 
N(J=l )^10xN(J=0), as o bserved in H2-detected QSO- 
DLAs llLedoux et al. 2003^ and in the Magellanic Clouds 

ijTumlinson et al 20021 . i.e. taking log N(J=0)=14.5 and 

log N(J=1)=15.5 (which is actually the above case 1), 
f should be less than or of the order of 10"^. Although 
our spectra have a lower spectral resolution than those 
normally used to study DLA systems along QSO lines of 
sight, the large H i column density in GRB 030323 allowed 
us to obtain an upper limit which is similar to the limits 
found in QSO-DLAs. We also examined the GRB 030323 
spectra for presence of absorption lines fror a vibrationally 
excite d molecular hydrogen predicted by IPraine fc Had 
1)2002!) . but these are also not de t ected. 

As shown by iLedoux et al.l l)2003|) . the lack of H2 
molecules in DLA systems is mainly due to the low metal- 
licity of the gas in addition to its particular physical 
conditions (density, temperature, UV flux). In particu- 
lar, H2 is usually not detected whenever the metallic- 
ity [X/H]<~1. In GRB 030323, the sulphur metaUicity, 
[S/H]=-1.26±0.20, is low enough to explain the lack of 
H2. An alternative explanation is that H2 close to the 
GRB has been dissociated by the strong UV/X-ray emis- 
sion; however such an emission would also ionize a large 
frac tion of the gas with w hich H2 molecules are associated 
(see iDraine fc Haoll2002|) . 

A dust depletion factor (i.e. the abundance difference 
between a dust-depleted element such as Fe or Cr and a 
non-depleted element such as Zn or S) of 0.2 dex at a 
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^ Observationally, [S/Si]=0 in gaseous absorbers when there 
is neither dust depletion nor ionization effects. 



2 3 
Redshift 

Fig. 6. Comparison of t he metallicities o f a sam ple of 
QSO-DLAs, taken from IProchaska et al.l l)2003a|) (open 
circles), with the two GRBs for which a metallicity 
has been determined (solid triangles): GRB 000926 and 
GRB 030323 (this paper). The GRB hosts are located at 
the metal-rich end of the QSO-DLA distribution. 
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metallicity of -1.26 (cf. [Fe/H]=-1.47±0.11) is also con- 
sistent with measu rements in QSO-DLAs (see Fig. 12 of 
iLedoux et al]l2003|) . However, this result is different from 
the an alysis of three GRB host galaxies bv lSavaglio et al.1 
l|2003(l . who find that the GRB host dust depletion is much 
larger than it is in QSO-DLAs. 

In Fig. we compare the metalliciti es (from Zn, S or 
Sil of a sample of QSO-DLAs taken from lProchaska et al.l 
" 2nn8a(l . with the GRB-DLAs for which a metallicity 



has been determined: GRB 000926 and GRB 030323 (this 
paper). For G RB 000926, we have adopted the value 
[Zn/H]=-0.25 of ICastro et al.l (2003), which is consistent 
with the curve-of-grqwth analysis value - [Zn/H]=-0.13 
- of ISavaglio et al.l l)2003|) . Since the neutral hydrogen 
colu mn density determin ation for GRB 000926 is not se- 
cure llFvnbo et alJl2nnil) . we assume an error of 0.3 dex. 
Although only two GRBs have measured metallicities. Fig. 
El suggests that G RB host galax i es are more metal rich 
than QSO-DLAs. ISavaglio et all l|2003(l already pointed 
out large Zn column densities in three GRB host galaxies 
(for which only GRB 000926 has a measured H i column 
density) with respect to QSO-DLAs, while they found the 
Fe column densities to be similar to those of QSO-DLAs. 
Hence, [Zn/Fe], a measure of the amount of dust depletion, 
is very large in their sample of GRB hosts with respect to 
QSO-DLAs, suggestive of a high dust content. Although 
we do not have an estimate of the Zn column density, 
the quantity [S/Fe] is a sim ilar measure. For GRB 03 0323, 
find [S/Fe]=0.2, while IProchaska eraP l|2003b(l find 



we 



<[S/Fe]>=0.4 (based on three systems in their sample 
for which this quantity is not an upper or lower limit). 

Using the measured metallicity and Hi column den- 
sity, we can check the low optical extinction that 
we inferred from the aft erglow photometry. Following 
IProchaska fc Wol^ l)2002l) . we assume that Av—Rvk 
N(Hi)host/4.9xl02i, where Ry =Av/E(B-V) is the total- 
to-selective extinction; Rv(MW)=3.1, Rv(LMC)=3.2, 
and Rv(SMC)=2.9 (see |Pe|[l992). The dust-to-gas ra- 
tio, K=10[^/"](1-10P''/^1), corresponds to the dust-to- 
gas ratio of the dust responsible for the extinction. 
The value (4.9±0.3)xl0^^ cm~^ mag~ ^ corresponds to 
the G alactic value for N(Hi)/E(B-V) ^Diplas fc Savagel 
Il994|) : for the LMC and SMC, we assum e the val- 

jKoornneei 



cm 

-2 



mag 

l|Bouchet et al. 



1985D, 



1982) and 



ues (2.0±0.5)xl0^ 
(4.4±0.7)xl02^ cm"^ mag 
spectively. With [S/H]=-1.26 and [Fe/S]=-0.21, we find 
K=0.02, and Av(MW)=0.08 mag, Av(LMC)=0.02 mag, 
and Av(SMC)=0.01 mag. These values are all consistent 
with the upper limits derived from the afterglow photom- 
etry (see Sect. I^J. 

As can be seen in Fig. 01 fine-structure lines of both 
Cii* and Sill* are detected. The Sin* lines have never 
been clearly detected in QSO-DLAs, which suggests that 
their origin is associated with GRB 030323, or that they 
can be found only in regions with very high neutral gas 
densities. Sill* A1264 has also been observed along the 
GRB 010222 sightline (I. Salamanca, private communica- 
tion) . The population of the fine-structure levels is a func- 



tion of the density of the absorbing medium and the am- 
bient photon-flux i ntensity ( Bahcal l fc Wolj Q,968). Using 
the calculations of ISilva fc Viegasi l|2002l) . we can make 
a rough estimate of the Hi volume density using the 
two un-saturated Sill* lines AA1309,1533 (the Cii* and 
Sill* A1264 lines are saturated), for which we measure log 
N(Siii*)~14.5. Assuming that [Siii/H] = [Sii/H], i.e. Si is 
undepleted onto dust grains, we obtain log N(Si 11)^^16.2 
(as log N(Sii)=15.84 and the Solar abundance difference 
between S and Si is 0.34) and n3/2/ni/2 '~ -1.7. This ra- 
tio corresponds to a volume density of uhi ^ 100 cm^'^ 
when the free electron density (uo) is assumed to be 10% 
of that of the Hi density (uhi), and uht ^ 10^ cm~ '^ 
when Uo < 10^^ uhi (see Fig. 8 of lSilva fc ViegaFJl2nn2^ . 
If these fine-structure lines originate in the same region 
as the neutral hydrogen, then the Si ii* medium is mostly 
neutral and the free electrons will mainly come from ion- 
ization of neutral atoms with an ionization potential lower 
than 13.6 eV, whose s olar abundance relatiy e to hydrogen 
is typically 10~ (see ISilva fc VieJ^l2002ll . This would 
result in an expected ratio Ue ^ 10""* nni. We have as- 
sumed that the fine-structure levels are populated by colli- 
sions between particles, and not through direct excitation 
by infra-red photons (although this mechanism is prob- 
ably not important i n the case of Sin*), or fluorescence 
I Silva fc Viegasll2'002l) . Under this assumption, we can di- 
vide the column density by the volume density to obtain 
an order of magnitude estimate of the size (diameter) of 
the absorbing region: ^ 5pc (taking uhi = lO'^ cm~^). As 
a comparison. Galactic molecular cloud sizes range from 
roughly 0.5pc t o SOp c l|Solomon et al.l 1798711 . Following 
ISilva, fc Vieg'^ l|2002|l . we can also estimate the mass of 
the Sin* absorbing cloud, M = mpN(Hi)lgjjj^ (where mp 
is the proton mass and Isi n* is the diameter of the Si ll* 
absorbing region) to be M = 2 x lO'' M©. However, the 
size and mass estimates would be upper limits if the Si ii* 
ions are only partly associated with the entire H i column 
and/or Sin column. For instance, if half of the Si ii absorp- 
tion would not be related to the Sin* absorbing region, 
the actual volume density would be roughly twice as large, 
the corresponding cloud size twice as small and the mass a 
factor of four smaller. If, on the other hand, fluorescence 
plays a non-negligible role, the size and mass estimates 
above would be lower limits. 

7. Lya in emission 

As can be seen in Fig. in the Lya trough Lya 
in emission is detected. Using splot within IRAF, we 
measured the center of the line at 5316.9±0.7A, cor- 
responding to z=3.3736±0.0006, with a FWHM of 237 
km s~^. The velocity offset with respect to the metal- 
absorption lines is -|-151±46 km s^^, i.e. the Lya emit- 
ting region is red-shifted with respect to the material 
responsible for the absorption lines. We measured a 
Lya flux of F = (1.2±0.1) x lO^^"^ erg cm^^ g-i in the 
1400V spectrum. The emission line is also detected, al- 
beit barely, in the lower resolution 300V spectrum with 
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Fig. 7. HST/ACS F606W imaging of the field of GRB 030323. On the left panel the 1' x 1' field surrounding the burst 
is shown, whose location is indicated by the cross hair. On the right panel the central 5"x5" region is shown, where 
in this case the image has been median filtered using a box with 3 pixels on a side. The position of the early afterglow 
and its uncertainty (0'.'06 in radius) is depicted by the black circle. At a distance of 0'.'14 (Ikpc), an extended object 
is detected. This is the probable host galaxy of GRB 030323 for which we measure AB(F606W)=28.0±0.3 mag. 



F = (1.0±0.3)xl0"i'^ erg cm^^ g-i^ Assuming Ho = 
70kms-iMpc-\ flu = 0.3 and f^A = 0.7, z=3.372 cor- 
responds to a luminosity distance of d^— 9x10^* cm. This 
transforms the observed flux into a Lya luminosity of l^hya 
= (1.2±0.1)xl042 erg s'^. 

Adopting a relation between measured Lya luminos- 
ity and star-formation rate of L/SFR = I C^^ erg s~^ per 
1 M0 (jKennicutt IQol lCowie fc HJ1998|) . the Lya star- 
formation rate in the GRB 030323 host galaxy is roughly 
1.2 Mq yr~^. This value has not been corrected for ex- 
tinction and hence is a lower limit. The Lya luminos- 
ity is roughly twice as large as the one obtained for the 
host of GRB 971214 a t z=3.42: (0.66±0.07)xl042 ej.g g-i 
ijKulkarni et alJll99^ but is at the lo w end of a saniple o f 
10 Lya emitters found at z ~ 3.4 bv lCowie fc R\l l|l998l) 
in and around the HDF and Hawaii deep field SSA 22 
(with LLya ranging from 1.2 to 8x10^^ er g s~^). However, 
a deeper survey by iFvnbo et alJ ll2003bl) has resulted in 
42 confirmed Lya emitting galaxies at z ~ 3, with most 
luminosities ranging from 0.4-2x10''^ erg s^^. 

Lya emission has also been o bserved in the troughs of 
half a dozen of QSO-DL As (e.g. |M 0ller fc Warren 'l99^ 
Pettini et all lUgS.: .Diorgovski et aL .1996: M&Wex et al.| 
1998HFvnbo et al.ll999l:lMcjller et alJ2002b|) . It is believed 
that this emission originates in the DLA host galaxy itself, 
and not in QSO photo-ionized regions when Zabs ~ zqso- 
In the GRB 030323 case, it is clear that the emission is 
produced by photo-ionization by massive stars (not neces- 
sarily related to the GRB) in the host galaxy. An origin in 



the immediate environment of the GRB is not possible, as 
this emission would be absorbed equally well as the after- 
glow continuum emission around Lya by the high neutral 
hydrogen column density along this sightline. 

8. HST imaging of the host galaxy 

The field of GRB 030323 was observed for 4 x 480s with 
HST/ACS in the F606W filter on July 20, 2003, starting 
at 23:00 UT. The dithered exposures were drizzled with 
the multidrizzle routine* to produce an output image with 
a scale of 0'.'033 per pixel. Fig. □ shows the I'xl' field 
of GRB 030323 on the left panel, and on the right panel 
the central 5"x5" region. The close-up image has been 
convolved with a median filter of 3 by 3 pixels, excluding 
the central pixel of the kernel in the median calculation. 
The position of the early afterglow has been projected onto 
the HST image using the 3min F0RS2 V-band image of 
March 26.35 (see Table|21); 8 objects were used to perform 
the transformation, with a resulting positional accuracy 
of 0'.'06. This error circle is shown on the zoomed image. 

At 4.3 pixels, or 0'.'14, an extended object is de- 
tected, which can be identified as the probable host of 
GRB 030323. The Gaussian FWHM of the point spread 
function of this object (5.0 pixels) is significantly greater 
than that of stars in the field (2.8 pixels), and SExtractor 
((Bertin fc Arnout^ll996l) classifies it as a galaxy, with the 

see ,http: / /stsdas. stsci.edu/pydrizzle / multidrizzle| 
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star-galaxy classification flag equal to 0.03 (which nor- 
mally is for a definite galaxy and 1 for a definite star). 
Adopting Ho=70 km s'^ Mpc~\ f}M=0.3 and f}A=0.7, 
0'.'14 at z = 3.372 c orresponds to an angular diameter 
distance of 1 kpc (see lHogdll999l) . 

For the host galaxy, we measured a magnitude 
AB(F606W)=28 .15±0.11 mag using SEx tractor's isopho- 
tal flux estimate ( Bertin fc Arnoutgll996|) , while we obtain 
AB(F606W)=28.1±0.2 mag with aperture photometry us- 
ing an aperture radius of 5 pixels. With a 10-pixel aper- 
ture radius, the magnitude increases to AB(F606W)=27.8 
mag. We adopted a zeropoint of 26.51 mag for the con- 
version of counts to AB magnitudes. The error estimates 
given above only include the Poisson errors. Assuming 
that the host galaxy has a flat spectrum in AB (i.e. /3=0 
in F^ oc ly^), the AB magnitude that we finally adopt: 
AB(F606W)=28.0±0.3 mag, is the same when converting 
it to the Johnson V band. 

From the UV continuum emission of the host galaxy we 
can obtain another crude estimate of the star-formation 
rate, independent of the SFR inferred from Lya in 
emission (see Sect. 0). The magnitude AB(F606W)=28 
mag corresponds to a flux Fj,=2.3xl0^'^^ erg s~^ cm~^ 
Hz~^, which at the luminosity distance of GRB 030323, 
di=9xl0^^ cm, results in Li,=5.4x 10^'' erg s_'^ Hz"^ 
Us: 



the SFR-L^(1500-2800A) relation of iKenniciittI 
Ijl998j) . this transforms to a SFR(UV) of 0.7 M© yr'^. 
We note that we have neglected Lya forest absorption at 
the blue end of the F606W filter (4700-7200A, while Lya 
is at 5317A), which would increase the SFR estimate by 
roughly 25% to around 0.9 Mq yr~^, which is very similar 
to the estimate from the Lya emission line (1.2 Mq yr^^). 



9. Discussion 

Our observations show that GRB 030323 occurred behind 
a very high H i column density, in an environment (imme- 
diate and host-galaxy combined) having a low molecular 
hydrogen fraction (f^lO~^), a low metallicity ([S/H]=- 
1.26±0.20) and a low dust content (k=0.02 ). For the DLA 
host of GRB 020124, iHiorth et all l|2003a^ also find evi- 
dence for a large H i column density with a low reddening. 
The inferred low dust content may be interpreted as a 
selection bias: GRBs that would occur in a dusty host 
galaxy would be harder to detect because they would 
then be fainter. However, in this case one would expect 
to observe many GRB afterglows with considerable ex- 
tinctioi i in the optical, for wh i ch there is no clear evi- 
dence l|Galama fc WiierslbOOll: iLazzati et al.ll2002a) . In 
apparent contradiction with this is the detec tion of sev- 
eral hos t galaxies in t he radio and sub-m m (e.g. 'Frail et al.' 
L2 QQ2; Barnard et a ll 12003: Berger et al. 2003). suggest- 
ing that at least some GRB hosts are dusty, as expected 
when most of the star forma tion in the univer se occurs 
in submm-bright gala xies fsee iRam irez-Ruiz et aL. 2002.). 
Dust destruction fe.g.lFruchter et alJl200lHDraine fc Had 



l20n2HPerna et aljl200 3!). which has been proposed to ac- 
count for the apparent discrepancy between the low op- 



tical extinction and lar ge (X-ray) gas column densities 
llGalama fc Wiiersll200lh . could play a role, but is not re- 
quired by our data: the reduced metallicity and hence the 
low dust-to-gas ratio in the host of GRB 030323 is suf- 
ficient to explain the combination of a larg e Hi column 
density with a low optical extinction (see also lHiorth et all 
l2003ali . 

Lya in emission is detected for GRB 030323, and we 
inferred a star-formation rate of about 1 Mq yr~^, which 
is in good agreement with the SFR value that we ob- 
tained from the UV c ontinuum emission of the host galaxy. 
iFvnbo et all l|2003a|) note that Lya is commonly observed 
in all GRB host galaxies at high redshift for which it could 
be detected. In contrast, only 25% of the Lyman-break 
galaxies are Lya emitters wi t h an equivalent width EW 
> 20 A llShanlev et a,lJl200,'^. iFvnbo et a,lJ ll2003ah sug- 
gest that this difference is due to GRB hosts having a low 
metallicity and a low dust content, consistent with our ob- 
servations of GRB 30323 and with those of GRB 020124 
ijHiorth et al.ll2003al) . We note that QSO-DLAs also have 
a low metallicity and a low dust content, but they rarely 
show Lya in emission. However, since most galaxy coun- 
terparts of QSO-DLAs are very faint, Lya in emission 
is not expected to be dete cted in most case s with the 
current detection limits fsee lFviibo et alJll999l). The low- 
dust inference for GRB 020124 ((Hiorth et al.ll2003al) and 
GRB 030323 is different from the results of ISavaglio et all 
(2003) , who find evidence for a high dust content in three 
GRB host galaxies. 

From the fine-structure lines Sin* AA1309,1533 , we 
estimated the Hi volume density of the gas producing 
this absorption: uhi = 10^ — 10^ cm^^, under the as- 
sumption that these fine-structure levels arc populated by 
collisions, and not through direct excitation by infra-red 
photons (which is not an important excita tion mechanism 
in th e case of Sill*), or fiuorescence fsee ISilva fc Viegii^ 
12002!) . This volume dens ity is higher than in ferred for 
QSO-DLA environments ijSilva fc Viegasll2002l). but typ- 
ical of Galactic molecular cl ouds (e.g. iBlitz fc Williamj 
I1999I: iReichart fc Pried l2002() . As these fines has never 
been clearly detected up to now in QSO-DLAs, the de- 
tection of these Sill* lines in the GRB 030323 spectrum 
suggests an origin in the vicinity of the GRB place of 
birth (e.g. the star- forming region in which it exploded). 
Combining the measured H i column density with the or- 
der of magnitude estimate of the H i volume density, we 
obtain a size (diameter) of ~ 5pc (taking uhi = lO'^ cm^'') 
and a mass of ~ 2x10'^ M0 for the Sin* absorbing region. 

With the volume density so high, one would ex- 
pect hydrogen molecules to be present, which, surpris- 
ingly, we do not detect. We obtain a rather strong 
upper limit on the mean molecular fraction of the 
gas in the GRB environment and the host galaxy: 
f=2iV(H2)/(2Ar(H2)+A^(Hi))<10-6. This could be ex- 
plain ed by the low metallicity of the gas fsee lLedoux et all 
l2003l) . but it may also be that the molecules in the GRB 
environment have been dissociated by the strong GRB 
UV/X-ray emission (e.g. iDraine fc HaQii2002.1 . In the lat- 
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ter case, however, the UV/X-ray flash wo uld also ionize 
the n eutral gas in the GRB vicinity (see iDraine fc Had 
which would make the high H i column density de- 
tection improbable. Therefore, a large fraction of the Hi 
column density may not be located close to the GRB ex- 
plosion site, but elsewhere in the host galaxy, while the 
high volume density Si ii* region (and the expected molec- 
ular hydrogen), is located in the vicinity of the burst. In 
this case, the disks of GRB host galaxies need to be much 
denser than the Galactic disk, as 7 random sight lines 
through the disk toward the location of the Earth would 
not resu lt in 5 H i column densiti es above 10^^ cm~^ (see 
Fig. 5 of lDickev &: Lockmanll99o|) . as is observed for GRB 
sightlines (see Fig. 2J. Finally, the population of the Sin* 
levels may have been partly caused by fluorescence of pho- 
tons from the GRB itself, in which case the volume density 
estimate above is a strict upper limit. 

The H I column density that we inferred toward 
GRB 030323 is higher than that of any (QSO- or GRB- 
) DLA measured using Lya in absorption. It is generally 
assumed that the apparent H i column density limit of 
N(Hi)~1022 atoms cm'^ for QSO-DLAs is due to an ob- 
servational bias against the detection of such high-column 
density systems, as these would obsc ure the background 
QSO if they contain s ome dust (e.g. lOstriker fc Heisleil 
Il984l iMLkPslEiii)- Ho wever, a radio-select ed QSO 
survey for DLA systems by lEllison et al.l l)200l[l did not 
uncover a previously unrecognized population of N(Hi)> 
10^^ cm~^ DLA systems in front of faint QSOs . An alter- 
native scenario was proposed bv ISchav3 l(200ll) : the lack 
of high Hi column density systems could be due to the 
conversion of Hi to H2 as the neutral gas density in- 
creases. This picture is consi stent with observ ations of 
Galactic molecular clouds fe.g. lBlitz Sz Williamail999.') . In 
GRB 030323, however, we do not find any evidence for the 
presence of H2 in addition to Hi to support this scenario. 
Future GRBs with possibly even larger H i column densi- 
ties than that toward GRB 030323 could provide further 
constraints to the existence of a rapid conversion of H i to 
H2 at high Hi column densities. 

We compared the metallicities and Hi column densi- 
ties of the (still very small) sample of GRB-DLAs with 
QSO-DLAs, and we found both quantities to be higher in 
GRB-DLAs than in QSO-DLAs. This is not surprising, as 
GRBs are now k nown to probe massive-star forming re- 
gions l|Stanek et al., 2003. iHiorth et aLll2003b^ where the 
gas density and the metallicity are higher than along ran- 
dom QSO sight lines through foreground galaxies. A KS 
test applied to the column densities shows that the prob- 
ability that the GRB- and QSO-DLA samples are drawn 
from the same parent distribution is very low (0.0006). 
On the other hand, two GRB afterglows have very low 
column densities. A large sample of high-resolution spec- 
tra of GRB afterglows could provide statistical informa- 
tion about the distribution of the gas in high-redshift star- 
forming regions, in addition to the evolution of the metal- 
licity and dust and H2 contents of GRB host galaxies. Such 
a sample can be created in the years to come thanks to 



rapid and accurate GRB localizations from future satellite 
missions such as Swift^ and EXIST^. 
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